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PACS 61 .30. Hn - Surface phenomena: alignment, anchoring, anchoring transitions, surface- 
induced layering, surface-induced ordering, wetting, prewetting transitions, and 
wetting transitions 

PACS 61.30.Dk - Continuum models and theories of liquid crystal structure 
PACS 64 . 70 . M — Transitions in liquid crystals 

Abstract. - Motivated by recent experimental observations [U. Delabre et al, Langmuir 24, 3998, 
2008] we reconsider an instability of ultrathin nematic films, spread on liquid substrates. Within a 
continuum elastic theory of liquid crystals, in the harmonic approximation, we find an analytical 
expressions for the critical thickness as well as for the critical wavenumber, characterizing the 
onset of instability towards the stripe phase. Comparing theoretical predictions with experimental 
observations, we establish the utility of surface-like term such as an azimuthal anchoring. 
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Introduction. — The formation of spatially periodic 
patterns in thin nematic films, spread on liquid substrates, 
is one of the well-known phenomena in physics of liq- 
uid crystals [IH21- Nematic liquid crystals (LC) are de- 
scribed by the director n, a unit vector, characterising 
an averaged preferred orientation of the molecules, which 
can vary throughout the sample [Hj- The liquid substrate 
tends to align this director n in the plane of the sample 
(planar anchoring), whereas the air interface imposes an- 
other preferred orientation, the one perpendicular to the 
interface (homcotropic anchoring). In thick nematic films 
(h > 1/im), the director n reorients in the vertical plane 
of the sample in order to relax antagonistic boundary con- 
ditions. However, in nematic films with submicron thick- 
ness, the competing interfaces start to "feel each other" 
resulting in undulations of the director. The in-plane sym- 
metry breaking is manifested in the formation of periodic 
patterns, like stripes, chevrons, zig-zags, squares [3]|4l[6][7] . 
An example of stripes in 6CB nematic film spread on wa- 
ter is shown in Fig. [T] The wavelength of the observed 
stripes 2 < L < 200 /im is typically two orders of mag- 
nitude larger compared to the thickness of nematic film 
20 nm < h < 0.5 /urn The presence of the free in- 

terface and small ratio of h/L oc 0.01 distinguishes our 
problem from the extensively studied Freedericksz transi- 
tions, where liquid crystals form periodic domains in pres- 
ence of electric and magnetic fields (SJUDJ- Moreover, our 




Fig. 1: Polarised-light microscopy image 1.8 x 1.3 mm of pe- 
riodically distorted 6CB/water LC at 27°C. The thickness of 
nematic film is h w 0.2 fim and the periodicity of stripes, the 
distance between two stripes, is L ~ 40 ^m> h. The alterna- 
tion of colour (dark and light) is associated with violent reori- 
entation (~ 7r/2) of the director n in the plane of the image. 



system is unique, since the wavelength L of stripes is con- 
trolled by the thickness h of nematic film in a non-trivial 
way. Therefore, one may think of the thickness "playing 
the role of the magnetic/electric field" in the Freedericksz 
transitions. 

First experimental observations of long-wavelength pe- 
riodic director distortions in thin nematic films were re- 
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ported in 1990 by Lavrentovich and Pergameshchik [T| 
for 5CB/glycerol LC system. Later, the authors pro- 
posed a model [21 E], applied to describe experimen- 
tal data of stripe domains with characteristic wavelength 
5 < L < 200 /xm as a function of the film thickness 0.14 < 
h < 0.49 /jm. A reasonable agreement between theory and 
experiment is achieved only by accounting for the splay- 
bend surface- like elastic term, with the value of clastic con- 
stant A"i 3 = -0.205An [3], where K n = 6.2 TO" 12 N [12]. 
Nevertheless, the predicted value for the lower threshold 
h c i = 0.14 /xm is almost one order of magnitude higher 
than h c i = 20 nm reported in 0] for 5CB/glycerol and 
in HI EH for different LC systems. Although the present 
studies [3" 1| 4" [| 6"] agree on the persistence of the striped phase 
up to the thickness of h C 2 « 0.5 /im (upper threshold), 
value for the lower threshold, h&, as well as the mecha- 
nism governing the instability remains an open question. 
The proposed theoretical models [2Tf4l |TTllT3"] suggest sev- 
eral mechanisms for the formation of stripes, known as 
A22, K24 or A24-A13 mechanisms. These mechanisms 
play a role only when the anchoring energy is weak and 
does not exceed the value of 10~ 6 J/m 2 as was suggested 
in [TT 1") . If we know the elastic constants and the anchor- 
ing energies by experimental measurements then we can 
choose the most suitable model to predict the formation 
stripes. However, the problem was addressed mostly nu- 
merically, assuming a priori either idealised values of elas- 
tic constants, e.g. A 24 = —K22 [2] and A 2 2 = An [4"IIT3"]. 
or the hybrid aligned cell (HAN) as only possible ground 
state [SlfTTj , which does not allow to understand the phys- 
ical picture as a whole. 

Recent experiments with different LC systems, 
such as 5CB/glycerol, 6CB/glycerol, 6CB/water, 
MBBA/glycerol [5J, suggest that stripes occur for films as 
thin as h c i ~ 20 nm, with the characteristic wavelength 
L ~ 2 /im. To quantify these observations we aim 
at reconsidering the problem by using an analytical 
framework, different from the one developed in [IH5] and 
identifying the critical thickness as well as the critical 
wavenumber, which determine the onset of the instability 
from a planar state towards a stripe state. Our paper 
concerns the systematic study of static instabilities within 
the Frank-Oseen continuum model of liquid crystals, 
assuming A 33 = Kn [E] together with A13 = [IB] , 
and supplemented by polar and azimuthal anchoring 
energies [5[. Dynamic instabilities have been previously 
studied in a more restricted context in [T7lfT8] . Our model 
is a first step for future generalisations. 

The influence of the azimuthal anchoring on the thresh- 
old thickness was studied in the pioneering work of Spar- 
avigna et al [5] . The authors have shown for the first time 
that decreasing the azimuthal anchoring and / or the value 
of the twist elastic constant A22 favours the formation of 
stripes. However, in that paper the role of the saddle- 
splay elastic constant A24 was disregarded, and only the 
case AT 2 4 = —K22 was considered. Taking into account 
both the azimuthal anchoring and the saddle-splay elastic 



constant would allow to compare theoretical results with 
experimental findings and to understand the picture as a 
whole. In the present paper we show that azimuthal an- 
choring, even being vanishingly small, restricts the undula- 
tions of the director and breaks the in-plane symmetry. As 
a result, at the critical point we find a well-defined wave- 
length for the stripes, which otherwise is infinite. Identify- 
ing a finite critical wavenumber and the critical thickness 
at the instability threshold makes possible a substantial 
comparison between theory and experimental data [6j[7]. 
Let us present first the model. 

Theoretical framework. — The nematic liquid crys- 
tals are described by an elastic free energy, quadratic in 
the director derivatives, given by [5] 

T eX = iyV{A n (V,n) 2 + A 22 (n,V x n) 2 + 
A 33 |nxVxn| 2 - (A 24 + A 22 )V[n(V, n)-(n, V)n]}, (1) 

where An, A22, A 33 and A24 are elastic moduli, char- 
acterising splay, twist, bend, and saddle-splay, respec- 
tively. Being a divergence, the last saddle-splay term can 
be transformed to a surface integral. Although being ne- 
glected for thick samples or samples with a strong anchor- 
ing boundary conditions, the surface-like terms are crucial 
in understanding the cause of instabilities in thin nematic 
films [SlimniHn] • The Erickscn inequalities [2T], given by 

K u > 0, A 22 + A 24 < 2 A n , |A 24 | < A 22 , (2) 

establish the relationships between the coefficients of the 
free energy ([1}, which guarantee the stability of the uni- 
form ground state with the director n = const. The 
anchoring energy is usually described by the Rapini- 
Papoular potential [22] and is written in terms of polar 
angle 9 and azimuthal angle (p as 

T a = ^ sin 2 ^ - 6 X ) + ^f- sin 2 (# 2 - h)+ 

+ ^ sin 2fl 2sin 2 ( ^_ v 5 2)j (3) 

where Wei and W02 are the polar anchoring strengths 
on the liquid substrate and at the air interface, re- 
spectively, W V 2 is the azimuthal anchoring strength at 
the air interface, with Wei > W02 W V 2- More- 
over, on the liquid substrate the preferred direction is 
9 1 = ir/2 (planar anchoring), whereas at the air in- 
terface £>2 = (homcotropic anchoring). Without loss 
of generality, the coordinate system is chosen in such 
a way that the preferred in-plane orientation at the air 
interface (p2 = 0. In the cases considered experimen- 
tally j3[l4f6] . i.e. 5CB/glycerol, 6CB/glycerol, 6CB/watcr, 
MBBA/glycerol, we are dealing with a weak anchoring 
regime, namely Wg oc 10~ 5 J/m 2 , the anchoring extrap- 
olation length being Lei = Kn/Wei — 0.35 /im and 
Lq2 = A11/W02 — 0.7 /tm. The azimuthal anchoring is 
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usually neglected at both interfaces, because the interfaces 
are isotropic and homogeneous. However, we believe that 
in the case of strong undulations of the director orienta- 
tion (see Fig.[l} considered here, azimuthal anchoring may 
become essential. A non-zero contribution to the free en- 
ergy originates from the tendency of n-CB molecules with 
long carbon tails to align along a common direction, thus 
penalising the perturbation of the director. 

In the following, we use the two constant approxima- 
tion Kn = -K33 = K. Then two ground states, yielding 
the minimum of T e i + T a , are the undistorted planar state 
with n = (1,0,0) or = tt/2 and distorted hybrid aligned 
state (HAN) with n = (sin0(z), 0, cos8(z)), is a polar 
angle varying along the thickness of the film h as 0(z) = 
01 + (02 — 0i)z/h. The normalised free energy density is 
given by W (0i, 2 ) = (0i -0 2 ) 2 + h cos 2 0i/Li+/isin 2 2 /L 2 . 
If the quadratic form w(0i, 2 ) is positive definite the pla- 
nar state (0i = 2 = 7r/2) is stable with respect to the 
HAN (0i ^ 2 ) state. The anchoring transition, between 
these two states happens at the Barbcro-Barbcri critical 
thickness [33] 

h c = L g2 -L gi , (4) 

which is approximately 0.35 fim for 5CB/glycerol LC sys- 
tem. However, if we do not require the director n to stay 
in the xz-planc, but allow a configuration corresponding to 
the minimum of the free energy, other equilibrium struc- 
tures might occur in response to the change of thickness, 
adopted (and not imposed) by the system. In the fol- 
lowing sections we explore a possibility of the transition 
between a homogeneous planar state and a stripe phase 
with periodic distortions of the director n in y-direction 
(see Fig. IB). 

Variational problem. — To analyse a relative stabil- 
ity of the planar ground state we consider the first varia- 
tion of the free energy. The perturbation of the director 
n = (1, 0, 0) can be written as n' = n + <m = ( sin(7r/2 + 
tf>) cos </>, sin(7r/2 + ip) sin </>, cos(7r/2 + ?/>)) , where <f> and ip 
are the small variations of the azimuthal and polar angles, 
respectively. The planar ground state is stable if the dif- 
ference in the free energy AJ 7 = ^{n + i5n} — J-{n} > 0, 
V <5n. If there exists n' critical for AT 7 , so that AT 7 < 0, 
then a distorted state is an equilibrium one Assum- 
ing d/dx = 0, since we are interested only in periodic 
modulation along the y-direction, and working in the har- 
monic approximation, we find from (JTJ), ([3]) the variation 
of the free energy as sum of the bulk fs and the surface 
fs contributions 



AJ- 



K 



(5) 



2L J dy { f3 + l dzfB }> 
f B = (#) a + tii'yf + {4>' y f + - 2r^ z , (6) 

<A 2 (0) i, 2 {h) , <t?(h) , 



fs = 



Lei L$2 L V 2 

-2(l- P ){m<Py(0)-^Wy(h)), (7) 



fi- z - t x 2 f + r X g' z = 0, 
WL ~ X 2 9 ~ rxfs = 0, 



where t = K 22 /K 7 r = 1 - t, p= (K 22 + K 24 )/K, L v2 = 
K/W v2 and L is the period of the stripes in y-direction. 

In equilibrium, the first variation of the free energy (|3"jl 
vanishes, 5(AJ 7 ) = 0. Thus, the functions ip(y, z) 
and <j){y,z), extremising AJ 7 , can be found from the 
Eulcr-Lagrange equations. Looking for a periodic so- 
lution written as ip{y, z) = f(z) sin(qy) and <p(y, z) = 
g{z) cos(qy) [TT], the Euler-Lagrangc equations for / and 
g are 

where the dimcnsionless variables z = z/h and X = qh 
are introduced. A general solution of these two coupled 
linear differential equations of second order (Eqs. ©) can 
be written in the form 

f(z) = A\ sinh(xz) + A 2 cosh(xz)+ 

+ A 3 S s'mh( X z) + A4S cosh( X z) 7 ([9^) 

g(z) = Bi sinh(x^) + B 2 cosh(x5)+ 

+ B 3 z sinh(x^) + B^z cosh(xz), ([Sb) 

with coefficients Ai , Bi dependent on four integration con- 
stants Ci, C 2 , C3, C4, given by the following recursive 
relationships 



Ai 



(l+t)Ci+rxC 4 



r(C 2 
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(l + t)C 2 -T X C 3 



r(d 



2x 



2t 



A 2 

A 4 
B 2 

Ba 



C, 



-(Ci + xQ) 
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Note that for t = 1 (t = 0), system ([5]) decouples, yielding 
the coefficients A3 = A4 = B3 = B4 = 0. 

Since, in the linear approximation, the total free en- 
ergy ([5]) is a quadratic form in the distortions <f> and tp, 
according to ([9]), (|T0]) it is also quadratic in the integra- 
tion constants Ci. Substituting the solution of the Euler- 
Lagrangc equations in ([5]), and integrating over z and y 
coordinates, we can rewrite AJ- in the matrix form as }20) 



A^ = \ ]T U.jCCj. i,j = 1, 2, 3, 4. 



(11) 



The coefficients My of the matrix M. can be also defined in 
a simpler way, using integration by parts over z as follows 



<9AF ^ „ r _ 



g/g 

dip 



d g/g 

dz dip' z 



dip 



dfs_ 

dip 



dfe 



z=0 



da 



Euler-Lagrange 

dfs , g/g 
dip 



z—h 



dipjh) 
da ' 



lower boundary 



upper boundary 



(12) 
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The variable tp is used for / and g, in order to simplify the 
form of UU). Since the functions f(z) and g(z) both satisfy 
the Euler-Lagrange equations ©, the expression under 
the integral in (|T2|) is zero. The coefficients Ci in (jTTj) can 
be found by extremising AJ- with respect to C^, namely 
dAF/dCi = MijCj = 0, with the latter equation being 
equivalent to the linear combination of the boundary con- 
ditions [20j (two boundary conditions for each function). 
A non-trivial solution of the system J^j ^ij Cj = exists 
if and only if the determinant of the matrix M. vanishes, 



(a) 



h, pm 



det M = 0. 



(13) 



In the case detA4 > 0, the planar state with C, = 
is stable, otherwise, an instability towards a periodically 
deformed (striped) state occurs. Within a harmonic ap- 
proximation there is no need in finding the coefficients Ci , 
because the instability threshold is completely determined 
by the matrix M. . In the next section we find the solution 
of the governing equation det Ai = 0, which defines an 
implicit relationship between the thickness h of the film 
and the wavenumber x, given that other physical param- 
eters of the LC system are known. The minimum of the 
curves allows to identify both the critical thickness and 
the critical wavenumber. 

The lower threshold. — We first consider an insta- 
bility of the planar state towards a stripe phase, neglect- 
ing the azimuthal anchoring contribution (W v 2 = 0) to 
the surface free energy 0. By substituting the solution 
of the Euler-Lagrange equations © into © , ■ an d us- 
ing (|12p. we find the coefficients M^- and consequently the 
determinant of the matrix A4 as 



det M = 



(I-t)V-(l+0 2 sinh 2 X 



{4/i 2 i 2 (sinhx) 2 



L0iLg 2 t 2 X 2 

+ h{Lei - L e2 )pt X [2p(l - t) X - (p - (4 - p)t) sinh(2 X )] + 

£nWY[p 2 (l-t) 2 X 2 - (P- (4-^) 2 (sinh X ) 2 ]}. 

(14) 

Figure ®l shows a stable solution of the equation 
det M(h,x) = for a typical value of t = 0.6 [51 H2] 
and Lqi — 0.35 /im, Lq2 = 0.7 /im [7J[T5], correspond- 
ing to 5CB on glycerol at room temperature. The mini- 
mum value of h on these curves yields the critical thick- 
ness h c \ , when the non-zero perturbations of the director 
n may appear. The critical wavenumber Xci at /i = h c i 
turns out to be zero irrespective of the value of p, mean- 
ing that the perturbations have an infinite wavelength. 
This finding contradicts experimental observations 4 , 6 
of a finite wavelength L « 2 /im at the lower thresh- 
old thickness. Nevertheless, slightly above the threshold 
(h > h c i) the model predicts a finite and small dimension- 
less wavenumber x- The expansion of (|14p around Xci = 
gives detM = -16h{h - h cl )t x 2 /(L ei L 62 ) + 0(x 4 ) with 
the critical thickness 

hd = (L e2 - L 01 )(2 - p)p, < p < 2* < 2, (15) 
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Fig. 2: The solution of the equation detTVf = for the free en- 
ergy (a) without azimuthal anchoring (|14|) . (b) with azimuthal 
anchoring as in (|16l) . The typical values Lg\ — 0.35 /im, 
Lg2 = 0.7 fim and t — 0.6 are chosen. The curves correspond 
to different values of p = t + K24/K, satisfying Ericksen in- 
equalities @, and 77 = L V 2/Lo2- (a) The minimum value of 
h(x) defines the critical thickness h c i (|15|) achieved at x = 0, 
V p. (b) Two symmetric minima with non-zero Xa and h& 
separated by a maximum h c i — Lg2 — Lei = 0.35 /im at X — 
for L V 2 > L" 2 — 0.455 pm (77 > 0.65). In degenerate case of 
vanishingly small azimuthal anchoring {L^2 — > 00) the critical 
parameters for (a) and (b) at p = 0.05 coincide. 



where p = t + K24/K satisfies the Ericksen inequalities 
Below h c i, the homogeneous planar state is stable 
(detM > 0), whereas above h c i the periodically mod- 
ulated stripe phase is an equilibrium one (det A'! < 0). 
Let us point out that i) the solution (|15[) gives the min- 
imum of /i(x)j since the second derivative d 2 h/dx 2 <x 
p(p — 2t)/(Lex — L02) + 0(x 2 ) > 0, is positive and ii) 
another formal solution for the critical thickness, h c i = 0, 
is stable only for p > 2t, which violates the Ericksen in- 
equalities. In the experiments on different LC systems [6], 
h c x = 40±20 ran, which is theoretically achievable only for 
a small values of p, when K24 ~ — K22 — — 0.6K. This fits 
in the range of the values —0.6 AT < K24 < 0.6 A, measured 
experimentally for 5CB LC compound [23]. However, in 
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(a) 

Xci 



(b) 




Fig. 3: Critical threshold h c \ and Xci = 2nh c i/L as function 
of p — t + K24/K < It calculated for different r\ = L V 2/L$2- 
The curves correspond to the local minima of det A4 (|18p . The 
values for t = 0.6, Lei = 0.35 /im, Lq2 = 0.7 /im are chosen. 
The long-dashed line in (b) corresponds to the case of vanishing 
azimuthal anchoring [Lnp2 — ► 00) and is given by (JXSJ . Strong 
azimuthal anchoring (77 oc 1) suppresses the formation of stripes 
Xci =0 (a), as predicted in [2]. 



the case of an exact equality, namely i\24 = (p = 

or p = 2t), the second derivative d 2 h/dx 2 vanishes, re- 
sulting in a saddle-point rather than a minimum for the 
curves h(x)- Notice that, according to (JT5J) , the critical 
thickness does not depend on the value of the twist elas- 
tic constant t, as was mentioned by Pergameshchik as 
well. The plot of (fT5jl , shown in Fig. [3j is in agreement 
with calculations of Sparavigna et al in [3]. In conclu- 
sion, the considered model provides a finite value for the 
critical thickness, however, at the critical point x is zero 
(infinite wavelength), which contradicts experimental ob- 
servations [6] . Our result is robust and cannot be improved 
by taking into account the if 13 surface-term, which would 
modify (JT2J) and would result in the rescaling of the critical 
thickness, achieved again at x — [25]. 

Now we aim at identifying the critical thickness and the 
critical wavenumber, taking into account the azimuthal 
anchoring W v2 on the nematic-air interface. In essence, 
the expected symmetry breaking should appear in the 



form of the determinant of the matrix A4, given by 

detM= (1 -^ 2 - ( r 1 + f 3 Sinh2X - 
L 81 L 92 L V 2t z X 

[h 3 t[(l + t) sinh(2 X )-2(l-t)x]+ 

+ h 2 x[4,t 2 (L gi + ^(sinhx) 2 - L e2 (coshx) 2 )+ 
+ L ei p(p(l - tfx 2 - (1 + t)(P - (4 - p)O(sinhx) 2 )] + 

+ hpt X 2 [2(L gi {Le 2 + L v2 ) - L e2 L ip2 )p(l - t)\+ 
+ (Lei(L e2 - L v2 ) + L e2 L v2 ) (p - (4 - p)t) sinh(2 X )] + 

L gi Lg 2 L v2 p 2 X 3 [p 2 (l-t) 2 X 2 - (P- (4-p)i) 2 (sinhx) 2 ] }. 

(16) 

Indeed, det M. = becomes a cubic equation with respect 
to h, rather than a quadratic one when W v2 = 0. The so- 
lution of this equation is plotted in Fig. [5b for different val- 
ues of the azimuthal anchoring W v2 ~ 3-10~ 5 -^lCP 7 J/m 2 
(L v a ~ 0.21 ~ 70 fan or 77 = L v2 /L e2 = 0.3 -r- 100), 
which is an unknown parameter for the studied LC sys- 
tems 07] . Depending on the value of rj the curves exhibit 
either a single minimum at {/i c i,Xci} = {Lg 2 — L#i,0} 
or two symmetric minima at x — iXci 7^ 0. In terms 
of the Landau theory of a second-order phase transitions, 
one might consider the wavenumber as an order parame- 
ter, being zero in homogeneous planar state and non-zero 
(or small) in a stripe phase, and the azimuthal anchor- 
ing playing the role of control parameter (e.g. temper- 
ature) The solution at Xci = looses its stability 
when d 2 h/dx 2 \h~>Le2-Lei < or 

L V 2 > L c ; 2 = — -i— ^- • [L 2 2 (l - 2t)+ 

v 3(L ei - L e2 )(l - p)H 

L 2 gi {l + 1 - 3p(l - P + t))+ L ei L g2 (t -2 + 3p(l - t))] . 

(17) 

In the case p = 0.05, illustrated in Fig. [2b, L™ 2 ~ 
0.455 /im. To study the behaviour of the curves around 
the critical point we expand (|16[) in powers of x as 

det M ^ a a {h) + a 2 {h) X 2 +a i {h) X i + 0{ X Q ), (18) 

where a = -h 2 (h + Lei — Lg 2 )/(L ei Lg 2 L ip2 ), the coeffi- 
cients a 2 and a.4 contain bulky expressions, which is not 
necessary to write explicitely. Unless 04 is always positive, 
the next order term should be added to the series (|18p. 
When a 2 > (equivalent to L v2 < L" 2 ([17])). the equilib- 
rium solution is Xci = 0, and h c i = Lg 2 — Lei is given by 
the condition det M. — ao — 0. For a 2 < 0, two symmetric 
minima occur at Xci = (2 fl 4), yielding the follow- 

ing condition det M — a — a2/(4a 4 ) = for the critical 
thickness h c i- The solution of these equations is shown 
in Fig. [3] assuming the values Lgi and t for 5CB/glycerol. 
The behaviour of the curves in Figs. [2J and [3] for other 
LC systems should be qualitatively the same, because it 
results from an intrinsic symmetry of the problem and it 
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is not caused by a particular choice of parameters. The 
plotted equilibrium values for Xci and h c \ correspond to 
the local minima of detyVf (fI7o]) . shown in Fig. [21 The em- 
ployed expansion of det M. is justified only if x J$ 1, which 
holds for small azimuthal anchoring in the whole range of 
p (see Fig. [3]). Note that in the limit of the vanishing az- 
imuthal anchoring (L V 2 — ► co) the solution for the critical 
thickness h c \ converges to (JTSJ) . The qualitative as well 
as the quantitative agreement between theory and exper- 
iment is achieved for p oc (i^24 oc —K22) and small, but 
non-zero, value of W V 2 oz 1CU 7 J/m 2 . 

Concluding remarks. — In this paper we considered 
the onset of instability of planar ncmatic films, subjected 
to competing boundary conditions, towards a periodically 
deformed state (stripes). In the harmonic approximation, 
we found an exact solution for the variational problem 
and were able to identify the critical thickness as well as 
the critical wavenumber, characterising the lower thresh- 
old instability. The analysis is performed within the con- 
tinuum theory framework, which does not necessarily hold 
for films with thickness of the tens of nanometers. Nev- 
ertheless, by taking into account the azimuthal anchoring 
at the nematic-air interface, we found a reasonable agree- 
ment between theory and experiment [Sl[7|. Considering 
only the saddle-splay surface term, turns out to be in- 
sufficient to identify a finite wavelength for stripes at the 
lower threshold. One could study the described phenom- 
ena above the threshold and decide to analyse the exper- 
imental data as in [3]- However, this would require a so- 
phisticated non-linear analysis for the free energy. Finding 
the upper threshold can give additional insight into exper- 
imental data and will allow to discuss the possibility for 
the formation of stripes in thick ncmatic samples observed 
in [T5] . The study of the instability from a hybrid aligned 
ncmatic towards stripe phase will be fulfilled in the future 
within the same framework. 
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